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ABSTRACT 

We perform several high resolutfon N-Body/SPH simulatfons of fow surface bright- 
ness galaxies (LSBs) embedded in cold dark matter halos to study how likely is bar 
formation in such systems. The behavior of various coUisionless galaxy models is stud- 
ied both in isolation and in the presence of a large perturbing satellite. We also consider 
models with a dominant gaseous component in the disk. We find that in general bar 
formation requires disk masses at least a factor of 2 higher than those inferred for LSBs 
under the assumption of a normal stellar mass-to-light ratio. Low surface density stel- 
lar disks contributing less than 10% of the total virial mass are stable within NFW 
halos spanning a range of concentrations. However, a purely gaseous disk can form a 
bar even for quite low masses and for realistic temperatures provided that cooling is 
very efficient (we adopt an isothermal equation of state) and that the halo has a very 
low concentration, c < 5. The bars that form in these LSB models are shorter than the 
typical halo scale radius - their overall angular momentum content might be too low to 
affect significantly the inner dark halo structure. Once formed, all the bars evolve into 
bulge-like structures in a few Gyr and can excite spiral patterns in the surrounding 
disk component. The recently discovered red LSBs show significant non-axisymmetric 
structure and bulge-like components and share many of their structural properties 
with the final states of our LSB models with massive disks. Our results imply that a 
bulge-like component must be present in any low surface brightness galaxy that ever 
went bar unstable in the past. 

Key virords: galaxies: dynamics — galaxies: evolution — galaxies: halos — methods: 
N-Body simulations — cosmology: theory 



1 INTRODUCTION 

Deep surveys carried out in the last decade or so have un- 
veiled a vast population of low surface brightness disk galax- 
ies (LSBs) that range from dwarf to fairly luminous galaxies 
and even giant spirals (Bothun, Impey & Malin 1986;Bothun 
et al. 1987; Schombert et al.l992; Impey et al. 1996). Typi- 
cal LSBs are late-type objects characterized by blue colors, 
low metallicity, very high gas fractions and low levels of star 
formation activity. More recently, a new class of gas-rich, 
late-type red LSB galaxies has been discovered (O'Neil et 
al. 2000) together with another population of red, early- 
type LSBs (Bejersbergen, de Blok & Van der Hulst 1999; 
Galaz et al. 2002). Some of the red LSBs have metallicities 
close to the solar value (Bergmann, Jorgensen & Hill 2002). 
Hence LSBs probably span a range of properties comparable 
to "normal" high surface brightness spiral galaxies (HSBs). 
The high gas fractions of LSBs suggest that they are quite 



unevolved objects with very low past and present star for- 
mation rates. The low metallicity (de Blok & Van der Hulst 
1998a,b) and the low gas surface densities (de Blok, Mc- 
Gaugh & Van der Hulst 1996a) can explain the low star for- 
mation rates of blue LSBs (Gerritsen & de Blok 1997; Van 
den Hoek et al. 2000). The fact that LSBs and HSBs seem 
to follow the same TuUy-Fisher relation (Zwaan et al. 1995; 
McGaugh et al. 2000) implies that the former have much 
higher dark matter contents than the latter within the re- 
gion typically probed by the observations (of order of a few 
disk scale lengths) . This conclusion stems from the low stel- 
lar mass-to-light ratios suggested by the blue colors of most 
LSBs. The latter imply that the low surface brightness is 
really the product of a low stellar surface density (de Blok, 
Van der Hulst & McGaugh 1996b) - a higher dark matter 
density is then required to explain why these galaxies have 
circular velocities comparable to those of similarly luminous 
HSBs (Verehijen & de Blok 1999). Hence LSB galaxies of- 
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for an ideal opportunity to study the distribution of dark 
matter and compare the observational results with current 
theories of structure formation. 

Predictions of cold dark matter simulations were found 
to disagree with the rotation curves measured from HI emis- 
sion in many dwarf galaxies and LSB galaxies (Moore 1994; 
Flores & Primack 1994)). Observed galaxies would have ha- 
los with a constant density core whereas simulated halos 
have cuspy density profiles falling as r^^ (the NFW model, 
Navarro, Frenk & White 1997) or even steeper (Moore et 
al. 1999a). While some observers have pointed out several 
limitations of the original measurements of rotation curves 
based on HI emission from atomic hydrogen and have thus 
analyzed several of the same galaxies using stellar Ha line 
(Swaters ct al. 2000; Van den Bosch & Swatcrs 2001; Swa- 
ters 2003), recent high resolution data obtained using both 
techniques seem to confirm that the majority of the rotation 
curves of LSB galaxies cannot be fit by the cuspy halo pro- 
files predicted by CDM cosmogonies (de Blok, McGaugh & 
Rubin 2001a,b;de Blok et al. 2001;de Blok & Bosma 2002). 
In dwarf galaxies strong supernovae feedback might play a 
role in affecting the overall mass distribution (Navarro et 
al. 1994), but the quiescent star formation histories of LSB 
galaxies together with the fairly large potential wells of many 
of them rule out such a scenario (Bell et al. 1994, 1999; 
Bergman, Jorgensen & Hill 2002). 

Recently Weinberg & Katz (2002) have suggested that 
the dark matter cusp might be erased thanks to the dynam- 
ical interaction with a stellar bar. It has been known for a 
while that a rotating barred potential would slow down due 
to dynamical friction against the halo background, shedding 
its angular momentum to the latter (Hernquist & Weinberg 
1994). Of course, the deceleration of the bar will be stronger 
in more massive halos. Debattista & Sellwood (1998, 2000) 
where the first to note that the dynamical interaction be- 
tween the bar and the halo could provide clues to the nature 
of the latter; they noted that if galEixies have massive dark 
halos as predicted by cold dark matter models it would be 
hard to explain why bars in many galaxies are quite fast ro- 
tators. On the other end, Valenzuela & Klypin (2002) claim 
that the net transfer of angular momentum to the bars is 
quite weak - bars slow down on a quite long timescale even in 
CDM halos as angular momentum is exchanged back from 
the halo to the individual stellar orbits supporting them 
once a sufficiently high force resolution is used in the simu- 
lations. Weinberg & Katz (2002) consider the idealized case 
of a very massive non-responsive bar and find that the res- 
onant transfer of energy and angular momentum from the 
bar to an NFW halo can actually change the density profile 
of the latter creating a constant density in only a few Gyr. 
The effect would occur only when the resonant orbits in the 
halo are well resolved, which requires a very large number of 
particles in an N-Body simulation, of order of a few millions. 
But Sellwood (2002), using a different numerical technique, 
finds the effect to be almost negligible for more realistic ini- 
tial setups of the galaxies in which the bar is not imposed 
from the start and can evolve instead of being just a rigid 
potential. Whereas further investigation of the effectiveness 
of the bar-halo interaction with self-consistent stellar bars is 
necessary, an even more basic question arises; can this mech- 
anism be at play in the most interesting case, namely that 
of LSB galaxies? How likely is that these galaxies would go 



bar unstable? Indeed, numerical studies on bar formation 
and bar-halo interactions have always employed models of 
high surface brightness galaxies. 

While bars seem to be ubiquitous among spiral galax- 
ies as a whole (Eskridge et al. 2000), low surface brightness 
galaxies are expected to be stable to bar formation due to 
a combination of low disk self-gravity and high dark matter 
contents (Mihos et al. 1997). But the current status of obser- 
vations seems to suggest a quite complex scenario. In fact, 
while blue LSB disk galaxies are indeed typically non-barred 
(although some dwarf LSBs have Magellanic bars), red LSB 
galaxies comprise several systems showing distortions, evi- 
dent bars and even bulge-like components (Bojersbergcn ct 
al. 1999). Tidal encounters might also trigger bar formation 
or other non-axisymmetric distortions even in LSBs (Mihos 
et al. 1997) and might actually be required to destabilize 
the gas in the disks and sustain star formation (Schombert 
ct al. 2001; Verde ct al. 2002). 

In this paper we study the formation of bars in models 
of LSB galaxies comprising a stellar and/or gaseous disk em- 
bedded in a dark matter halo whose structure is motivated 
by the results of CDM models. Such an approach is novel 
among studies of bar formation in galaxies, having been only 
partially adopted by Valenzuela & Klypin (2002). We have 
carried out a large set of high resolution N-Body/SPH sim- 
ulations using the parallel binary tree + SPH code PKD- 
GRAV/GASOLINE (Stadel, 2001; Stadel, Richardson & 
Wadsley 2002; Wadslcy, Stadel & Quinn, 2003). 

Weinberg & Katz (2001) argue that, although present- 
day low surface brightness galaxies might be stable to bar 
formation, they could have undergone a bar instability in 
the past when they first assembled a cold, massive disk. 
Although such a bar would have likely triggered a burst 
of star formation, this being difficult to reconcile with the 
many hints pointing towards a fairly smooth star forma- 
tion history for these galaxies (Bergman, Jorgensen & Hill 
2002), observations still caimot exclude that at least the red- 
der among these galaxies might hide a faded massive disk 
(O'Neil et al. (2000)). The bulge-like components seen in 
many red LSBs might be the result of secular evolution of 
an old bar (Combes 1991; CaroUo et al 2001). In brief, the 
origin and nature of LSB galaxies is still subject to debate 
and therefore we will consider a range of models covering a 
vast parameter space in terms of masses and internal struc- 
ture of disks and halos. Our only two constraints will be that 
the rotation curves of the models have to resemble those of 
observed LSBs, i.e. they have to be slowly rising out to sev- 
eral disk scale lengths (de Blok & McGaugh 1997), and that 
disks have scale lengths bigger than those of HSBs having 
the same luminosity (Zwaan et al. 1995). 

Models with a major gaseous disk can represent either 
LSBs during their early evolutionary stage or those numer- 
ous present-day LSB galaxies in which the baryonic mass is 
mostly contributed by the gas component, (de Blok & Mc- 
Gaugh 1997). We note that the same disk model realized 
as purely gaseous instead of purely stellar is not expected 
to have the same stability properties; it has been shown 
that fluid configurations tend to be more stable because 
pressure is generally isotropic while the corresponding stel- 
lar analog, the velocity dispersion, is generally anisotropic 
(Cazes & Tohline 2000). However, the study of the stabil- 
ity of purely gaseous disks embedded in dark matter ha- 
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los is new to our knowledge; works exist on the stability of 
generic axisymmetric or triaxial fluid configurations (Gazes 
& Tohline 2000, Barnes & Tohline 2001) and even on bar- 
unstable exponential gaseous galactic disks (Friedli & Benz 
(1993; 1995) but all these studios did not include any realis- 
tic dark matter component in their models and are therefore 
extremely idealized when applied to the dynamical evolution 
of real galaxies. The stability of gaseous and stellar systems, 
both uniformly rotating and differentially rotating, was also 
studied by Cliristodoulou ct al. (1995) - these authors tried 
to extend to gaseous disks within dark halos a criterion for 
stability against bar formation originally developed by Efs- 
tathiou. Lake & Ncgropontc (1982) for stellar disks but did 
not test their conclusions with numerical simulations. 

Finally, we will also consider the case of perturbations 
tidally induced by massive satellites, as expected during the 
hierarchical build-up of structures. Strong tidal interactions 
with even more massive galaxies would also occur in a hi- 
erarchical scenario and surely will drive bar formation but 
they would also transmute these fragile disk galaxies into 
spheroidals or even destroy them (Moore et al. 1999; Mayer 
et al. 2001a,b). 

The paper is organized as follows; in section 2 wo de- 
scribe the models and the initial setup of the simulations, in 
section 3 wo illustrate the results of the simulations, section 
4 contains the discussion and a summary follows in the last 
section. 



2 GALAXY MODELS AND SIMULATIONS 

Galaxy models are built as in Mayer et al. (2001b;2002) us- 
ing the technique originally developed by Hernquist (1993) 
(see also Springel & White 1999). We use a system of units 
such that G = 1, [M] = 6.5 x 10^ M© and [R] = 6 kpc. The 
models comprise a dark matter halo and an embedded stellar 
or gaseous disk (or both). Structural parameters wore chosen 
in order to obtain slowly rising rotation curves resembling 
those published for LSB galaxies (e.g. de Blok et al. 2001a,b, 
de Blok, McGaugh & Rubin 2001, de Blok & Bosma 2002, 
Van den Bosch & Swaters 2001). We start by choosing the 
value of the circular velocity of the halo at the virial radius, 
Vvir, which, for an assumed cosmology (hereafter flo — 0.3, 
A = 0.7, Ho = 65 kms~^Mpc~^) automatically determines 
the virial mass, Mvir, and virial radius, Rvir, of the halo 
(Mo, Mao & White 1998). We choose Vvir = 75 km/s. Halos 
have NFW density profiles (Navarro, Frenk & White 1995, 
1997) with different halo concentrations,c, and spin param- 
eters, A (Mo, Mao & White 1998). The concentration is de- 
fined as c = Rvir/vs, where is the halo scale radius; the 
spin parameter is defined as A = J\E\^^^G~^ M~^^'^ , where 
J and E are, respectively, the total angular momentum and 
total energy of the halo and G is the gravitational constant. 

The value of the concentration c basically defines what 
fraction of the total mass of the halo is contained within its 
inner regions, where the baryonic disk lies; the concentra- 
tion increases with decreasing mass and, for a given mass, 
has a scatter of roughly a factor of 2, mainly due to dif- 
ferent formation epochs (Bullock et al. 2001; Eke, Navarro 
& Steinmetz 2001). The average concentration for galaxies 
with Vvir ~ 75 km/s is ~ 12 in the standard LGDM model 
assumed here. However fitting rotation curves of LSBs with 



NFW profiles often requires c ^ 5 (Van den Bosch &: Swa- 
ters 2001; de Blok & Bosma 2002), at the lower end of the 
allowed range of values. We consider three values for the 
halo concentration, c = 4, c = 7 and c = 12. 

Placing galaxy disks inside cosmological halos would in 
principle require to have first solved the problem of galaxy 
formation in the LGDM scenario. Incidentally, whereas for 
years simulations of galaxy formation within CDM mod- 
els have been plagued by the so-called angular momentum 
"catastrophe" (Navarro & Steinmetz 2000), producing only 
tiny disks an order of magnitude smaller than those of real 
galaxies, new SPH simulations with considerably higher res- 
olution and an improved treatment of gasdynamics find this 
problem to be significantly alleviated (Governato et al. 2002, 
see also Sommer-Larsen et al. 2002; Thacker & Gouchman 
2001). These new simulations produce disks bearing scaling 
relations with their dark matter halos reasonably close to 
the predictions of the semi-analytical models by Mo, Mao & 
White (1998) - the latter are able to match the properties 
of spiral galaxies assuming angular momentum conservation 
of the baryons as they cool into the halos and settle into 
a centrifugally supported disk(Fall & Efstathiou 1980). We 
thus construct models of LSB galaxies that follow the scal- 
ing relations of Mo, Mao & White (1998). The procedure 
used to assign structural parameters is described in detail 
in Mayer et al. (2001b) and Mayer et al. (2002). Here we 
recall that we use exponential disks (with a few exceptions 
indicated below) and that their mass and scale length arc de- 
termined primarily by Myir and A and, to a minor extent, by 
the disk/halo mass ratio, fd, and halo concentration c (the 
latter two parameters contribute to specify the potential en- 
ergy of the system and thus the rotational energy needed for 
centrifugal support). The adiabatic contraction of the halo 
in response to the accumulation of baryons at the center 
is also taken into account (Springel & White 1999). As we 
mentioned in section 1, rotation curves of LSB disk galaxies 
suggest that these systems are extremely dark matter domi- 
nated, with stellar disk/halo ratios ~ 0.03 or less, namely at 
least a factor of 2 lower than those of HSBs of comparable 
luminosity (O'Neil et al. 2001, Ghung et al. 2002). However, 
this diff'erencc probably reflects the fact that in LSBs a larger 
fraction of the baryons, sometimes most of them, are found 
in the gaseous component; indeed the total disk mass (stars 
+ gas) relative to the halo mass is quite similar in LSBs 
and HSBs (McGaugh ct al. 2000). Previous works on the 
detailed mass modeling of LSB galaxies have indeed found 
a typical value of 0.065 for the total disk/halo mass ratio in 
these systems (Hernandez & Gilmore 1998). We thus adopt 
the view that LSBs are simply more extended than HSBs 
because of larger halo spin parameters and eventually lower 
halo concentration, as originally suggested by Dalcanton et 
al. (1997) and Mo, Mao & White (1998), (see also Jimenez 
et al. (1997) and Hernandez & Gilmore (1998)) and as con- 
firmed more recently by Jimenez, Verde & Oh (2003) in their 
modeling of Hq rotation curves. The disk mass fraction fd 
is usually set equal to 0.05 or less in our models but we also 
consider cases in which /d = 0.1, such a massive disk being 
expected in scenarios where LSBs have a massive faded stel- 
lar component (O'Neil et al. 2001;Clmng et al. 2002; Galaz et 
al. 2002). The highest value of fd is still consistent with the 
upper limit set for by nucleosynthesis, ~ 0.13 (Fukugita, 
Hogan & Peebles 1999); in this case, however, we are im- 
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plicitly assuming that nearly all the baryons ended up in the 
disk (it is somewhat an extreme assumption - see Verde et al. 
2002). Jimenez et al. (1998), who coupled a stellar popula- 
tion synthesis model to a simple galaxy formation scheme in 
the context of the hierarchical clustering, find that A > 0.05 
is required to match the surface brightness and colors of sev- 
eral "blue" LSBs. In our models, halos have spin parameters 
either 0.065 or 0.1, larger than the mean value found in cos- 
mological simulations, A ~ 0.04 (e.g. Gardner 2001, Lemson 
& Kaufmann 1998) - these values yield a disk scale length, 
Rh, in the range 2-6 kpc, consistent with observations (e.g. 
Zwaan et al. 1995, O'Neil et al. 1998). Although our halos 
have Vvir ~ 75 km/s, the different halo concentrations and 
the addition of the stellar disk produce a peak rotational 
velocity, Vpeak (at ~ 2Rh) between 95 and 130 km/s, within 
the range of the maximum rotation speeds measured for the 
majority of LSBs in the samples by de Blok & McGaugh 
(1997) and by de Blok & Bosma (2002). Indeed it is only 
Vpeak which is accessible to observations based on HI or Ha 
kinematics, not Vvir, which is too far out in the halo. The 
rotation curves of some galaxy models are shown in Figure 
1. 

The setup of the stellar disk is complete once even the 
Toomre parameter, Q{R), is defined (Toomre, 1964) This 
corresponds to fixing the local stellar velocity dispersion ctr, 
as Q{R) = aRK/3.36GT,s, where k is the local epicyclic fre- 
quency , G is the gravitational constant and Es is the disk 
surface density. The Q parameter at ~ 1 disk scale length 
will be given as a reference; the latter is typically set equal to 
1.2 (see Table 1 for the precise values adopted). Q, in com- 
bination with other parameters (Binney & Tremaine 1987; 
Mihos et al. 1997), determines the stability of the disks to 
bar formation, Numerical studies done in the past suggest 
that Q ^ 2 is needed for stability against stellar bar growth 
in isolated galaxies (Athanassoula & Sellwood 1986; Priedli 

2000) ; however these studies were conducted using galaxy 
models of "normal" HSB spiral galaxies, while our LSB mod- 
els have a lower disk self-gravity and could well be stable 
even for smaller Q values (Mihos et al. 1997). 

We note that models with c = 4 — 7 , A = 0.065 and 
fd ~ 0.05 have a B band central surface brightnesses, /iog, 
comprised between 22.5 and 23.5 mag arcsec^ for a B band 
stellar mass-to-light ratio (M/Lb)^ = 2. We avoid more ex- 
treme values of the surface density /brightness on purpose; 
indeed out first objective is to investigate if the bar instabil- 
ity is possible at least among LSB galaxies at the bright end 
of the surface brightness distribution (de Blok & McGaugh 
1997). The models with massive disks {fd = 0.1) would have 
a central surface brightness 0.5 mag lower than the thresh- 
old often used as a definition of LSB galaxies, /ios = 22.5 
mag arcsec"^, for a stellar mass-to- light ratio ~ 2, but, if 
representative of faded disks with M/L, ^ 5 (O'Neil et al. 

2001) , would have Hb > 23.5 mag arcsec"^. 

In the runs employing a gaseous disk this has a tem- 
perature of 7500 K; the kinematics of the neutral hydrogen 
in the disks of spiral galaxies yield typical velocity disper- 
sions consistent with this temperature (Martin & Kennicutt 
2001). The gaseous disk has either a constant or an expo- 
nential surface density profile with the same scale length of 
the disk (see Mayer et al. 2001b); in fact some LSB galaxies 
have rotation curves indicating that the HI component has a 
profile flatter than the stellar disk (e.g. de Blok & McGaugh. 
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Figure 1. Rotation curves of galaxy models with purely stellar 
disks out to three disk scale lengths. The solid line denotes the 
total curve, while the dot-dashed and the dashed lines represent 
the separate contribution of, respectively, dark matter and stars. 
The name of each model is indicated in the panels (see Table 1). 



1997). In a gaseous disk the Toomre parameter is defined as 
Q{R) = VaK/-nGT.g, where Vg is the sound speed and Eg 
is the surface density of the gas. For the assumed disk and 
halo masses a temperature of 7500 K implies Q ^1, namely 
comparable to that of the coUisionless runs. 



We do not include radiative cooling in the treatment of 
gasdynamics but we adopt an isothermal equation of state 
to model dissipation. This is somewhat idealized as it stands 
on the assumption that the thermal energy generated by, e.g. 
shocks and supernovae explosions, can be instantly radiated 
away, i.e. cooling is assumed to be very efficient, however 
simulations of galaxy formation and evolution which explic- 
itly include both radiative cooling and heating show that 
the temperature of the gas in galactic disks stays always 
close to 10** K (Gerritsen & Icke 1997; Katz, Hernquist & 
Weinberg 1992; Sommer-Larsen et al. 2002; Governato et al. 
2002; Hernquist & Katz 1989; Katz & Gunn 1991; Navarro 
& Steinmetz 2000; Governato et al. 2002 - note however 
that lower temperatures are not allowed in many of these 
works because they usually adopt cooling functions with a 
cut off at about 10* K). Supernovae might heat the gas 
very efficiently if the energy of their explosions is partly con- 
verted into turbulent motions (Thacker & Couchman 2001; 
Springel & Hernquist 2002) instead of being entirely con- 
verted into thermal energy, but their global impact on galax;- 
ies bigger than dwarfs is not yet established both observa^- 
tionally and theoretically (Martin 1999; MacLaw & Ferrara 
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Table 1. Parameters of the initial models. Column 1: Name of the model. Column 2: Disk/halo mass ratio. Column 3: halo concentration. 
Column 4: halo spin parameter. Column 5: disk gas fraction. Column 7: type of gas surface density profile. Column 8: £<£ stability 
parameter. Column 9: Toomre Q parameter at 1 disk scale length; Column 10: X2 parameter for swing amplification at 1 disk scale 
length. Column 11: mass of the satellite relative to total virial mass of the system. Column 12: whether the system forms or not a bar 
("tr" stands for transient) 
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fd 


c 


A 


fg 


gas profile 


£d 


Q 


X2 




bar 


Lmd0c4 


0.038 


4 


0.065 








1.123 


1.2 


2.6 





no 


Lmdlc4 


0.05 


4 


0.065 








0.945 


1.2 


1.9 





no 


Lmdlc4Q2 


0.05 


4 


0.065 








0.945 


1 
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no 


Lmdlc4Q3 
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4 
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0.945 


0.5 
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12 
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1.55 





yes 
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4 
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1 


cxp 


1 
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4 
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1 


exp 


0.703 


1.2 


1.36 





yes 


Lmdlcl2g 


0.05 


12 


0.065 


1 


exp 


1.02 


1.2 


2.47 





no 


Lmd2cl2g 


0.1 


12 


0.065 


1 


exp 


0.707 


1.2 


1.55 





yes 


Lmdlc4sg 


0.05 


4 


0.065 


0.5 


exp 


1.41 


2.4 


3.8 





no 


Lmdlc4gc 


0.05 


4 


0.065 


1 


const. 


1.4 


2.2 


10 





no 


Lmdlcl2gc 


0.05 


12 


0.065 


1 


const. 


1.2 


1.5 


13 





no 


Lmdlc4sat 


0.05 


4 


0.065 








0.945 


1.2 


2 


0.03 


no 


Lmdlc4Q3sat 


0.05 


4 


0.065 








0.945 


0.5 


2 


0.03 


tr 


Lmdlcl2sat 


0.05 


12 


0.065 








1.02 


1.2 


2.47 


0.03 


no 




5 10 15 5 10 15 20 25 

R[kpc] R[kpc] 




5 10 15 2 4 6 8 10 12 



R[kpc] R[kpc] 



Figure 2. Rotation curves of galaxy models with purely gaseous 
disks out to three disk scale lengths. The solid line denotes the 
total curve, while the dot-dashed and the dashed line represent 
the separate contribution of, respectively, dark matter and gas. 
The name of each model is indicated in the panels (see Table 1). 



1999; Benson et al. 2003). We recall that the general expec- 
tation is that LSB galaxies should be quite stable to non- 
axisymmetric instabilities; in this context the assumption 
of an isothermal equation of state will provide the most fa- 
vorable condition for the formation and survival of bars in 
gaseous disks by forcing the disk to remain cold (for instance 
it might underestimate heating when a disk has already en- 
tered a phase of strong instability) - in reality the same 
systems can only be more stable if they retain some of the 
heat generated during their dynamical evolution. 

However, in order to understand how sensitive are the 
results to the assumed equation of state, we evolve the same 
initial conditions with both an isothermal and and adia- 
batic equation of state (not shown in Table 1); the adia- 
batic runs represent the situation at the opposite extreme of 
the isothermal runs, namely radiative cooling is completely 
switched off. A final type of models comprise systems in 
which the stellar and gaseous component make an equal 
contribution to the disk mass - this can reflect an evolution- 
ary stage intermediate between those of models with purely 
gaseous and purely stellar disks. The rotation curves of some 
of the models with a gas component are shown in Figure 2. 

Gasdynamical simulations were carried out with GASO- 
LINE, an extension of PKDGRAV (Stadel 2001) which 
uses SPH to solve the hydrodynamical equations (Wads- 
ley, Stadel & Quinn, 2003; Stadel, Richardson & Wadsley 
2002). The gas is treated as an ideal gas with equation of 
state P = (7 — l)pu, where P is the pressure, p is the den- 
sity, u is the (specific) thermal energy, and 7 is the ratio of 
the specific heats (adiabatic index) and it is set equal to 5/3 
(we are assuming that the gaseous disk represents the atomic 
hydrogen component of the galaocy). In its general form the 
code solves an internal energy equation which includes an 
artificial viscosity term to model irreversible heating from 
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Figure 3. Color-coded face-on view of the stellar density after 
3 Gyr of evolution for some of the coUisionless galaxy models 
(see Table 1). Prom top to bottom; Lmdlc4, Lmdlc4Q3, Lmd2c4 
and Lmd2cl2 are shown. Brighter colors correspond to higher 
densities; a central bar is evident in some cases. Boxes are 20 kpc 
on a side 



shocks. The code adopts the standard Monaghan artificial 
viscosity as well as the Balsara criterion to reduce unwanted 
shear viscosity (Balsara 1995). In the isothermal case the 
thermal energy is constant over time and so no thermal en- 
ergy equation is required (any heating resulting from the 
artificial viscosity term, which is still present in the momen- 
tum equation, is radiated away). In the adiabatic case the 
thermal energy can rise as a result of compressional and 
shock heating and can drop following decompressions. 

We also explore the case of LSB galaxies tidally dis- 
turbed by a massive companion having a mass comparable 
to that of the LMC. Other regimes of external tidal pertur- 
bations should not be relevant here. In fact more massive 
satellites would rapidly merge at the center before being 
substantially stripped (Colpi, Mayer & Governato 1999; Taf- 
foni et al. 2003) and would likely destroy the fragile disks, 
while repeated fly-byes with even more massive galaxies in 
dense environments, like galaxy clusters or groups, would 
actually transform LSB galaxies into spheroidals or SO sys- 
tems (Moore et al. 1999; Mayer et al. 2001a,b). For this set 
of runs we use LSB galaxies with purely stellar disks and 
the satellite is a deformable spherical halo with a fairly con- 
centrated NFW profile (c=15), a high concentration making 
the satellite stiffer and the associated perturbation stronger. 
The models used for these runs are indicated in Table 1. 
The satellite moves in the same direction of the disk rota- 
tion (prograde), which can maximize the tidal effects due to 
resonances between its orbital frequency and the frequencies 
of disk stars (Velazquez & White 1997). We consider both 
a very eccentric (apo/peri=15) and a nearly circular orbit 
(apo/peri=2), both having a pericenter of 30 kpc, namely 
grazing the disk; indeed, while repeated tidal shocks can be 
effective in driving global instabilities (Mayer et al. 2001b), 
near-resonant interactions between the satellite orbital fre- 
quency and the internal frequency of the galaxy might also 
excite non-axisymmetric modes in the disk (Weinberg 2000; 
Weinberg & Katz 2002). 

For the models with a single-component disk we use 
500,000 or 1 million particles for the halo and 50,000 par- 
ticles for the disk (either gaseous or stellar). In the models 
with a two- components disk we use 25,000 particles for both 
components and the same number of halo particles as before. 
In order to test how our results depend on resolution, we also 
performed runs in which the latter is doubled or tripled in 
either the dark or the baryonic component. The resolution 
used for the halo is more than adequate to study physical bar 
formation - previous work showed that, with a few hundred 
thousand particles in the halo, spurious bar modes cannot be 
triggered in otherwise stable galaxies (Mayer et al. 2001b). 
The only other existing work in which N-body simulations of 
bar formation with similarly high resolution were performed 
is Valenzuela & Klypin (2002); compared to them we use a 
lower force resolution, i.e we employ a (spline) softening cor- 
responding to 0.05-Rh (for both the stars and the gas) while 
they use O.OlRh- However, a small softening is not necessar- 
ily a good choice - force resolution must be balanced with 
mass resolution in order to avoid increasing the noise due to 
the discrete representation of the systems (Hernquist, Hut 
& Makino 1993;Moore, Katz & Lake 1996) and our conser- 
vative choice has been widely tested in this respect (Mayer 
et al. 2001b). Yet, in order to test how are results depend on 
force resolution, we have also run selected models a second 
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Figure 4. EUipticity parameter of the stellar eomponent in some 
of the coUisionless simulations after 3 Gyr of evolution. The thick 
solid line refers to model Lmd2cl2, the thick dashed line to model 
Lmd2c4, the thin dashed line to model Lmdlc7, the thick long- 
dashed line to model Lmdlc4Q3, the thin solid line to model 
Lmdlcl2 and the thin long-dashed line to model Lmdlc4. 

time using the same softening as in Valenzuela & Klypin 
(2002) without finding any significant difi^erences in the re- 
sults. 



3 RESULTS 

We evolve the various initial galaxy configurations for up to 

10 Gyr, which corresponds to several disk dynamical times 
(this being of the order of a few times 10* years) and already 
represents a significant fraction of the cosmic time. The de- 
veloping of non-axisymmetric structure in our galaxy models 
can be seen in the projected color coded density maps shown 
in this section. We also quantify the strength of such distor- 
tions by means of the ellipticity parameter e = 1 — b/a, where 
a and b are the major and intermediate axis of the baryonic 
component of the galaxies (measured on the basis of the 
principal mass moments, e = if the galaxy is axisymmet- 
ric). The ellipticity has been found to be a good measure of 
the strength of bars in a large number of studies as it corre- 
lates with more physical parameters like the ratio between 
the radial and the axisymmetric components of the force 
(Laurikainen et al. 2002; Das et al. 2001; Laurikainen, Salo 
& Rautiaiucn 2002;Martin 1995;Regan & Elnicgrccn 1997). 
Although both bar mass and ellipticity would be needed to 
fully specify the bar contribution to the overall stellar po- 
tential of the galaxy, the latter is the only parameter easily 
accessible by the observations (Martin 1995). 

3.1 CoUisionless runs 

LSB models with either low or high concentration and 
fd < 0.1 do not undergo any bar instability if Q > 1 (taking 



Figure 5. Ellipticity parameter of the gaseous component in gas- 
dynamical runs. Models Lmd2cl2g (thick solid line for the isother- 
mal simulation and thick long-dashed for the adiabatic simula- 
tion), Lmdlc4g (thick dashed line for the isothermal simulation 
and thin dashed for the adiabatic simulation) and Lmd2c4sg (thin 
long-dashed line) are shown after, respectively, after, 3, 5 and 3 
Gyr (i.e. at the time where the strongest non-axisymmetry is ob- 
served in each of them). 



Q at roughly 1 disk scale length as a reference like in Ta- 
ble 1). Stability of these models should be even more robust 
at higher resolution as numerical noise is suppressed and 
spurious non-axisymmetric modes induced by the discrete 
representation of the coUisionless fluid gradually disappear 
(Hcrnquist 1993). A transient bar appears after 3-4 Gyr in 
LSB models with c — 4 halos and fd ^ 0.05 only when Q is 
as low as 0.5 (see Table 1). The bar scale length is small, 
roughly 0.5Rh- The results for such "light" disk models are 
shown in Figure 3 and 4. Instead, models with massive disks, 
fd = 0.1, become bar unstable regardless of halo concentra- 
tion and the value of Q. With such high disk masses bars 
are strong and long-lived (bottom panel of Figure 3, Figure 
4 and 10). Bars form rapidly, after about 2 Gyr, of order of 
a few disk dynamical times at i? = Rh, and are ~ IRh in 
low concentration halos, and much longer, ^ 2Rh, in high 
concentration halos. We tested that our results are not com- 
promised by numerical effects, especially two-body scatter- 
ing between stellar particles and more massive halo particles 
(Moore et al. 1996) and lack of force resolution that might 
bias the evolution of the galaxy potential (Dehnen 2001); 
we found good convergence when we compared a)runs with 
500,000 particles halos with runs having 1 million particles 
halos and b)when we compared the latter runs with runs 
having the same number of particles but a force resolution 
five times higher (i.e. a softening five times smaller) in both 
the dark matter and the stellar component. 

At higher values of the halo concentration a higher 
shear is present near the center and swing amplification 
of m = 2 modes is expected to be weak unless Q is ex- 
tremely low (Binney & Tremaine 1987). The susceptibility 
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of a differentially rotating thin disk to swing amplification 
of some particular mode m can be measured by the pa- 
rameter Xm = i?K^/27rmGE (Toomre, 1981). For m = 2 
modes, X2 < 3 is required for strong swing amplification 
in potentials associated with flat rotation curves (Binney 
& Tremaine 1987) but a lower threshold, X2 < 2, has been 
found for systems with slowly rising rotation curves as those 
of our LSB models (Mihos et al. 1997; Mayer et al. 2001b). 
The trend of increasing stability with increasing concentra- 
tion is consistent with the fact that the X2 parameter at 
R = Rh is nearly twice as big for halos with c = 12 com- 
pared to halos with c = 4 with having the same fd (see Table 

1) . X2 is indeed below the expected threshold for instability 
only for the models with fd = 0.1. Once initiated, swing 
amplification cannot proceed if a strong Lindblad inner res- 
onance is present that cuts off the feedback cycle (Binney 
& Tremaine 1987; Sellwood & Evans 2001). This can be an- 
other way by which galaxies can avoid to develop a bar. We 
checked the height of the Lindblad barrier for different mod- 
els. We conclude that stellar frequencies that can be swing 
amplified always lie above the inner Lindblad resonance and 
hence this mechanism is negligible in our simulations, al- 
though the barrier is certainly higher when a high concen- 
tration halo is present. In brief, we think that our results 
mostly refiect the structural difference between systems in 
which swing amplification can start and systems in which 
it cannot start and that such difference is due to a combi- 
nation of halo concentration and disk mass (both these two 
quantities enter the definition of X2 through, respectively, a 
and Ed). We observe a weak spiral structure even in mod- 
els with high concentration and fd ~ 0.05 and with up to 
a million particles (as shown by the non-zero ellipticity of 
the corresponding models in Figure 4). - however we cannot 
exclude that these features are partially driven by numerical 
noise (Hernquist 1993). 

We find that the £d parameter introduced by Efstathiou, 
Lake & Negroponte (1983) provides a good measure of the 
stability against bar format ion. This parameter is defined 
as Ed = Vpeak/^y {GM/Rh), where Vpeak (see section 2) is 
proportional to the mass of the dark halo within about one 
disk scale length, and the denominator is proportional to the 
disk self-gravity. In Table 1 we list the values of Sd for the 
galaxy models and we indicate whether the galaxy model 
is found stable or unstable to bar formation. In Efstathiou, 
Lake & Negroponte (1983) it was found that Sd 1.1 was 
required for stability against bar formation; we find that the 
criterion is somewhat weaker, Ed ^ 0.94. This slight discrep- 
ancy is likely due to the much lower resolution used in the 
N-Body simulations of those authors (more than an order 
of magnitude less particles than even the lowest resolution 
simulations presented here) which was likely enhancing bar 
formation. The fact that Ed does not account for several dy- 
namical aspects of bar formation, like the stabilizing effect 
of a large velocity dispersion of the disk and the possible in- 
terruption of swing amplification by the Lindblad barriers , 
and yet it provides a good fit to our results, strengthens our 
interpretation that the halo/disk mass ratio within the typ- 
ical disk radius, which is related to both halo concentration 
and disk mass fraction (see rotation curves in Figures 1 and 

2) is really determining the behavior of our galaxy models 
(had we included models with "hard centers" produced by. 
e.g., bulges, the incompleteness of this criterion would have 




Figure 6. Color-coded face-on view of the stellar density of model 
Lmdlcl2 (left) and Lmdlc4Q3 (right) perturbed by a satellite on 
a prograde orbit with, respectively, apo/peri = 2 and apo/peri = 
15 (see text). Brighter colors represent higher densities. Snapshots 
are taken after 6 Gyr, when the satellite has performed roughly 
three orbits. Boxes are 30 kpc on a side. 

been more manifest, see Sellwood & Wilkinson (2001)). Even 
so, Ed alone provides only a sufficient condition for stability 
to bar formation; a necessary and sufficient condition (at 
least for the present models) can only be obtained by cou- 
pling it with another parameter, for instance Q, and this is 
shown by comparison between models Lmdlc4, Lmdlc4Q2 
and Lmdlc4Q3 - for the same Ed a (transient) instability 
happens only when Q < 1 (see Table 1 and Figure 3). 



Galaxy models with highly concentrated halos are sta- 
ble even when satellites perturb them during several close 
passages over 10 Gyr (model Lmdlcl2sat, see Figure 6). 
Models with low concentration halos undergo a bar instabil- 
ity only if the disk has the lowest value of Q (Lmdlc4Q3sat) 
but this is weak and transient exactly as in the case where 
the galaxy is isolated, and the galaxy only shows mild non- 
axisymmetric distortions after several Gyr (Figure 6). We 
find these conclusions to hold regardless of whether the orbit 
of the satellite is nearly circular or nearly radial. On a very 
eccentric orbit the impulsive heating at pericenter drives the 
formation of outer spiral arms in the disk. Therefore, the 
perturbing satellite does not appear to trigger the bar in- 
stability in the disks: this is in contrast to what observed 
when an LSB satellite interacts with a much more massive 
galaxy halo (Mayer et al. 2001a, b), suggesting that only 
extremely strong perturbations can stimulate the growth of 
bar modes in these low surface density disks. 

In conclusion there seems to be only a small region of 
the parameter space possibly covered by LSB galaxies where 
the formation of a stellar bar is possible within a cuspy halo. 
In particular, these disks are stable within NEW halos as 
concentrated as typically expected in LCDM (c > 10). In 
order to grow bars, models need a baryonic fraction as high 
as 10 % of the total mass for a spin parameter ~ 0.065. This 
is more than a factor of 2 higher than the baryon fraction 
estimated for most observed LSB galaxies ((Hernandez et al. 
1998) and is higher by an even larger amount compared to 
the fraction of baryons estimated to be in stars for a normal 
stellar mass-to-light ratio (de Blok & McGaugh 1997). Of 
course one could imagine that current surveys are still miss- 
ing a fraction of the baryons of LSBs because these galax- 
ies are very extended and/or have a significant old stellar 
population not detectable in the optical (O'Neil et al. 2000; 
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Galaz et al. 2002); even so the required baryon content is 
quite high, being close to ^ 0.13, the upper hmit indicated 
by nucleosynthesis (Fukugita, Hogan & Peebles 1999 - but 
note that the new WMAP data allow for a sensibly higher 
upper limit, ~ 0.19 (Spergel et al. 2001)). Hence the mod- 
els with massive disks are somewhat extreme as they de- 
mand that nearly all of the baryons available to the galaxy 
have already cooled into the disk. As mentioned in section 
2, if (M/Lb)^ = 2 model Lmd2cl2g has a central surface 
brightness high enough to be close to the upper limit usu- 
ally adopted as a definition of an LSB galaxy in surveys, 
fxoB ~ 22.5 mag arcsec"^, while it would be a more typical 
LSB when placed in a low concentration halo ((/io_g '~ 23.5 
mag arcsec"^); this is because in our modeling more con- 
centrated halos naturally produce more compact disks (Mo, 
Mao & White 1998). One option to bring down the initial 
surface density/brightness (keeping Vvir and fd fixed) is to 
start from a larger value of the spin parameter, which would 
reduce the disk surface density by spreading the same disk 
mass over a larger radius. To explore the latter hypothe- 
sis we evolved a galaxy with A = 0.1 and c = 12 (model 
Lmd2cl2b, see Table 1) that would have a lower central 
surface brightness ( ^ 23.5 mag arcsec'^), near the average 
for LSBs. Such a model still goes bar unstable after ~ 2.5 
Gyr. Hence there is some flexibility in the choice of the disk 
surface density that can lead to bar formation, even for a 
highly concentrated halo, but always provided that the disk 
is sufficiently massive. 

3.2 Gasdynamical runs 

The simulations employing an exponential gaseous disks 
with fd — 0.05 and an isothermal equation of state show 
results that are sensibly different from those of the corre- 
sponding collisionless runs. A distinct bar forms within a 
c = 4 halo (Lmdlc4g, see Figure 7) after 5 Gyr, while only 
a short-lived bar-like distortion was observed in the cor- 
responding collisionless run and only for the lowest initial 
value of Q (model Lmdlc4Q3, see Table 1). Only milder non- 
axisymmetric patterns, either spiral-like or oval, appear after 
a few Gyr in halos having c = 12 (Lmdlcl2g). Later these 
patterns evolve into clump-like structures near the peaks of 
the surface density distortions. We note that the bar which 
forms in the Lmdlc4g is visibly shorter than that produced 
in any of the stellar-dynamical models (its radius is ~ O.SRh) 
and it also appears much later; we intend to explore in detail 
the formation path of this type of gaseous bar in a forthcom- 
ing paper. Once formed, the bar quickly evolves (in 1-2 Gyr, 
of order of its rotation period) into something resembling 
an oval distortion and then a bulge-like component (Figure 
8). It is thus a really short-lived feature during the evolu- 
tion of the gaseous disk. Rising the initial disk temperature 
by about 50% (model Lmdlc4gb in Table 1) is enough to 
prevent the growth of the bar. Nonetheless, this simulation 
suggests that cold gaseous disks appear more prone to un- 
dergo non-axisymmetric instabilities with respect to their 
stellar analog (see also Figure 5 on the ellipticity) . The rea- 
son is that the temperature is constant with time for model 
Lmdlc4g while its stellar dynamical analog, the velocity dis- 
persion, is not constant in the corresponding collisionless 
models (e.g. Lmdlc4), instead it actually increases with time 
, driving Q towards large values. For example, after 4 Gyr 




Figure 7. Color-coded face-on view of the stellar density of model 
Lmd2c4g. Brighter colors represent higher densities. Snapshots of 
both the isothermal run (top) and the adiabatic run (bottom) are 
shown at 5 Gyr (left) and 8 Gyr (right). Boxes are 30 kpc on a 
side. 



(just before the bar appears in model Lmdlc4g) the average 
velocity dispersion of model Lmdlc4 has increased by a fac- 
tor 2.5 within the disk scale length in response to weak spiral 
instabilities, and Q > 2 results throughout the disk, too high 
for the bar to develop. We tested whether the gaseous disk 
would respond in a similar way if cooling was inefficient; we 
thus evolved the same disk adiabatically for 6 Gyr and we 
found that, as expected, the bar does not form in this case. 
In fact, after 4 Gyr, the temperature has grown by more 
than a factor of 5 in the inner few kpcs (reaching 35000 K), 
which raises Q by a about a factor of 2 (Q oc T^/^). The 
disk now looks even smoother than the stellar disk, i.e. even 
less non-axisymmetric structure is present (compare Figure 
7 and Figure 3). However, the isothermal calculations are 
probably more realistic at least in a global sense as cooling 
should be efficient enough to keep the temperature of the 
HI disk below 10.000 K (Martin & Kennicutt 2001 - see also 
section 2). Therefore, we tend to conclude that in reality 
gaseous disks should be more susceptible to grow bars com- 
pared to stellar disks, although this does not imply anything 
on the relative longevity of the two types of bars (see be- 
low) . The disk evolution is also dependent on the type of the 
gas profile. In models where the gaseous disk has a constant 
gas-density profile (and is evolved isothermally) a strong bar 
does not appear neither in low nor in highly concentrated 
halos. In these models the surface density of the gaseous disk 
is too low everywhere (see the rotation curves in Figure 2 
and the values of the stability parameters in Table 1) and 
therefore it is not surprising that the growth of bar modes 
is suppressed. 
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Figure 8. Color-coded face-on view of the stellar density of model 
Lmd2cl2g. Brighter colors represent higher densities. Snapshots 
of both the isothermal run (top) and the adiabatic run (bottom) 
are shown at 3 Gyr (left) and 6 Gyr (right). Boxes are 30 kpc on 
a side. 



Massive gaseous disks {fd — 0.1, models Lmd2cl2g or 
Lmd2c4g) become strongly bar unstable in a similar fashion 
as their stellar counterparts (Figure 8). The bar instabil- 
ity appears even when the same models are evolved using 
an adiabatic equation of state. The instabihty in massive 
gaseous disks is "dynamical" as in the collisionless case - 
it arises early, on a timescale comparable to the dynami- 
cal time of the disk (Figure 8). The bar morphology evolves 
faster than in the stellar dynamical case, resembling more 
an oval distortion after just one rotation period (~ 3 Gyr, 
see Figure 8). This difficulty shown by the gaseous disks in 
sustaining a barred potential was also reported by Friedli & 
Benz (1993) who found that the periodic xi orbits that sup- 
port a stellar bar are quickly destabilized by shocks devel- 
oping in their loops. Shocks are indeed clearly visible along 
the bar edges in our gasdynamical runs and the elongation 
of the density distribution decreases rapidly, especially in 
the adiabatic runs, where a higher pressure develops that 
smears out the density perturbation (see Figure 5). 

Christodolou et al. (1995) have extended the stability 
criterion based on Sd to gaseous disks within halos, deriv- 
ing Ed ^ 0.9 as a condition for stability to bar formation in 
the latter systems. Their result was not based on numerical 
simulations of disk-f halo systems but was obtained from a 
comparison between quite idealized ellipsoidal fluid configu- 
rations and their stellar analogs. Our numerical results con- 
firm that gaseous models are slightly more stable than stellar 
disks if we compare models starting with a similar initial Q 
and Ed and there is no cooling (as in the adiabatic case). 
However, in general ~ 1 for stabihty, i.e. we obtain a cri- 
terion slightly stronger than for stellar disks if we take into 
account the effect of cooling (isothermal case). Nonetheless, 
the fact that gaseous bars appear to weaken sooner than 
stellar bars (this being also dependent on the equation of 




Figure 9. Color-coded face-on view of the stellar density of model 
Lmdlc4sg (see Table 1). Brighter colors represent higher densities. 
Snapshots are taken after 3 Gyr (top) and 7 Gyr (bottom). The 
stellar distribution is shown on the left, the gas on the right. Boxes 
are 30 kpc on a side. 

state) suggests that their overall dynamical impact on the 
galaxy (including the halo) might be less important. 



In the model where stars and gas contribute equally to 
the disk (model Lmdlc4gs)each of the two baryonic compo- 
nents has a surface density reduced by 50% compared to the 
corresponding models having a single baryonic component 
(Lmdlc4 or Lmdlc4g) and a bar does not form regardless 
of the halo concentration (Figure 9). Analytical calculations 
and numerical experiments have shown that a disk made of 
stars and gas can be more susceptible to gravitational insta- 
bility than either of its individual components would (Jog 
1992; Elmegreen 1995; Jog 1996); however, those results ap- 
ply when at least one of the components is marginally unsta- 
ble (yielding Q ~ 1), whereas in the present case the surface 
density is very low for both components and, as a conse- 
quence, all the relevant parameters have values well within 
the regime of stability (see Table 1). 

3.3 Morphological evolution of the bars 

In models that go bar unstable (both stellar and gaseous) we 
always observe a morphological evolution from a bar into a 
bulge-like central structure over a few bar dynamical times 
(a few Gyrs, see Figure 10 and 11), in agreement with the 
secular evolution scenario (Combes & Sanders 1981; Combes 
et al. 1990; Merritt & Sellwood 1994; CaroUo et al. 1999). 
In the latter vertical instabilities occur locally as a conse- 
quence of resonances between the bar rotation speed and 
the vertical oscillation frequency of the stars (Combes & 
Sanders 1981) or simply result because the orbital families 
supporting a fiat, radially anisotropic bar become unstable 
to bending modes (Pfenniger 1984; 1985; Raha et al. 1991; 
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Figure 10. Color-coded face-on view of tlie stellar density of 
model Lmd2cl2. Brighter colors represent higher densities. From 
the upper left to the bottom right, snapshots are shown at, re- 
spectively, 1,3,4.5,7 Gyr are shown. Boxes are 25 kpc on a side 




Figure 11. Color-coded edge-on view of the stellar density of 
model Lmd2cl2. Brighter colors represent higher densities. From 
the upper left to the bottom right, snapshots are shown at, re- 
spectively, 1,3,4.5,7 Gyr are shown. Boxes are 25 kpc on a side 



Sellwood & Merritt 1994). When viewed edge-on the galax- 
ies go from a disky shape to peanut-shaped soon after the 
bar forms and become progressively rounder (Figure 11). As 
the rotating bar evolves stars also lose angular momentum 
to the surroundings and lead to an increased surface density. 
The increase in central density is more pronounced in the 
gaseous disk (compare Figure 12 and Figure 13). In fact, in 
addition to the gravitational losses already present in col- 
lisionless systems, some dissipation of angular momentum 
can occur in the shocks near the bar edges; in the adiabatic 
runs this is compensated by the subsequent heating and ex- 
pansion of the gas and in the end the inner profile looks 
more similar to that of the collisionless runs. The bar per- 
sists for several Gyr, although it weakens as time goes on 
(the elongation of the density distribution diminishes, see 
last snapshot of Figure 10), maybe because of the growing 
central mass concentration (Hasan & Norman 1990; Friedli 
& Benz 1993; Norman, Sellwood & Hasan 1996). An early- 
type spiral results, but one in which the disk has still a very 
low surface density. The central bulge has an exponential 
profile with a scale length 5-10 times smaller than that of 
the surrounding disk (Figures 12, 13 and 14). The disk is still 
exponential but considerably flatter than at the beginning. 
Its scale length more than doubles, reaching as much as 9 
kpc in one case; the surface density decreases correspond- 
ingly and would lead to a disk central surface brightness a 
few magnitudes lower than in the initial conditions (Figure 
12, 13). These and other properties of the final states of our 
bar-unstable models closely match those of the red early- 
type LSB galaxies studied by Bejersbergen et al. (1999). In 
the next section we will discuss in detail how far the com- 
parison with such class of galaxies can be pushed. 



Our simulations do not include star formation. We 
might ask how long will a galaxy remain in a mostly gaseous 
state as assumed for some of our models. Kennicutt (1998) 
and Martin & Kennicutt (2002) have shown that the Toomre 
Q parameter gives a good indication of the local density 
threshold above which star formation can occur. Observa- 
tions of galaxies suggest that Q < 1.44 is required for star 
formation to proceed. Our gaseous disks satisfy the latter 
criterion nearly everywhere (except in the inner few hundred 
parsecs) but the star formation rates calculated from the lo- 
cal initial gas surface density as in Kennicutt (1998) would 
be as low as 0.3MQ/yr within 10 kpc (corresponding to 
a few disk scale lengths) for models with fd ~ 0.05. Such 
modest star formation rates fall near the lowest measured 
by Kennicutt (1998) for spiral galaxies and are comparable 
to the star formation rates that Gerritsen & de Blok (1999) 
obtain in their simulations of LSB galaxies. Therefore, such 
a galaxy will remain mostly gaseous for a good fraction of its 
lifetime. In addition, these low star formation rates imply a 
weak supernovae feedback with little impact on the thermal 
and dynamical evolution of the gas, this being consistent 
with our assumption of an isothermal evolution. Of course, 
when significant non-axisymmetric structure, bars or clump- 
like features form all these assumptions could break locally 
as the surface density grows above some threshold. In the 
light disk models {fd — 0.05) the surface density increases, 
on average, by only a factor 2 in the bar region (in the inner 
few kpcs in Figure 14). Instead, in models with heavy disks 
the average surface density within 1-2 scale lengths increases 
by a factor of 10 and would result in a typical star forma- 
tion rate of > 2OM0 /yr, comparable to that of actively star 
forming HSB galaxies (Bell et al. 1999). The highest star for- 
mation rates would occur in the center, inside the growing 
bulge-like component. As we mentioned, we believe that the 
final state of the heavy disks resembles the red LSB galaxies 
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studied by Bcjorsbcrgcn ot al. (1999). Clearly more quan- 
titative predictions ol the luminosity and color evolution of 
the simulated galaxies demand that we include star formar 
tion and the mechanisms that can regulate it, for instance 
heating by the UV radiation from hot stars and supcrnovao 
(Gerritsen & Icke 1997;Gerritsen & de Blok 1999). Simu- 
lations including these additional mechanisms will be the 
subject of a forthcoming paper. Nonetheless, our main fo- 
cus here is on the initial development of non-axisymmetric 
structure in both stellar and gaseous disks and in this re- 
spect LSB models start with surface densities low enough to 
justify neglecting star formation and its effects. 



4 DISCUSSION 

We have shown that bar formation in low surface brightness 

stellar disks is unlikely unless they are more massive than 
usually believed. Light disks generate only transient non- 
axisymmetric distortions even when starting from a state 
with Q < 1. Instead, purely gaseous disks can become bar 
unstable even for quite low meisses provided that their ha- 
los have low concentrations and their temperature stays 
low enough to maintain Q Si 1 (as shown in the simula- 
tions adopting an isothermal equation of state). Growing 
modes are damped more easily in stellar disks because the 
velocity dispersion increases in response, rising Q to values 
well above 1; gaseous disks behave similarly only when the 
cooling is completely switched off (like in the experiments 
with an adiabatic equation of state). Although adopting an 
isothermal equation of state might seem simplistic, the tem- 
peratures we need to keep Q ~ 1.2 over most of the radial 
extent of our models, T = TSOO-ft', is typical of the HI com- 
ponent of spiral galaxies (Martin & Kennicutt 2001). 



Concentrations for halos forming at a redshift > 2 in 
LCDM models are sufficiently low (c < 6 at 14 < 80 km/s) 
to include even the values required to drive bars in light 
gaseous disks (Bullock et al. 2001). The combination of a 
cold, mostly gaseous disk and a low concentration halo was 
probably not uncommon at high redshift and perhaps so 
were barred LSB galaxies. However, when we look at the 
present-day distribution of halo concentrations in cosmo- 
logical simulations (this is basically the combination of the 
distributions for objects formed at any redshift) we see that 
such low values are found for only a few percent of the sys- 
tems (Eke, Navarro & Steinmetz 2001; Bullock et al. 2001). 

Therefore an even stronger result of our work is that 
halos as concentrated as typically occurs in LCDM models 
(in passing we note that even the highest value of the halo 
concentration considered here, c = 12, is still conservative) 
inhibit bar formation, even when a satellite as massive as 
the Large Magellanic Cloud is perturbing the disk during 
several close passages. 

On the other end, massive disks only 10% lighter than 
the dark matter halo become bar unstable regardless of halo 
concentration but the resulting bar is longer (relative to the 
initial disk scale-length) and, possibly, stronger for higher 
halo concentrations (see section 3.1). This is consistent with 
the findings of Athanassoula (2002) and Athanassoula & 
Mistriotis (2002), who claim that a higher halo mass within 
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Figure 12. The final stellar surface density profile (thick solid 
line) of model Lmd2cl2 (see Table 1) is plotted against the intial 
purely exponential profile (thin solid line). Exponential fits to the 
inner (short-dashed line) and outer (long-dashed line) part of the 
profile are shown, with scale lengths of, respectively, 900 pc and 4 
kpc (while the initial disk scale length was ~ 3 kpc) . The surface 
density is expressed in the units of the simulation 



the disk region produces longer and stronger bars because 
these lose more efficiently angular momentum through reso- 
nant exchange with halo particles (but see also Athanassoula 
2003). 

As explained in section 3.3, bars evolve into bulge-like 
structures in a few Gyr in both stellar and gasdynamical 
runs. The bulge reaches a surface density significantly higher 
than the rest of the disk (Figures 12, 13, and 14); in models 
with gaseous bulges a mean star formation rate ~ IM^/yr 
would result in light disk models following Kennicutt (1998) 
- this implies fairly rapid gas consumption timescales, ^ 1 
Gyr. In contrast, the rest of the disk will continue to form 
stars with a rate from 10 to 20 times lower and will even- 
tually keep a sufficient gas reservoir to form stars until the 
present epoch, thus maintaining rather blue colors as those 
of many observed LSBs. Star formation rates between 10 
and 30 times higher would result in the bulges forming from 
massive disks; in this case a real starburst will likely convert 
most of the gas into stars in ~ IQ^years. In the latter case 
the bulge B band surface brightness would be comparable 
to that of bulges of normal early type spirals, ixb ~ 19 mag 
arcsec^, for a stellar- mass-to- light ratio ~ 1 (typical of a 
young stellar population) based on the stellar or gas surface 
density measured in our simulations. The bulge would even- 
tually fade following the burst; neglecting any further evo- 
lution of the surface density, if the stellar population under- 
goes passive evolution and reaches the typical stellar mass- 
to-light ratios of spheroids after a few Gyr, (M/Lb)^ ^ 3—4, 
the final central B band surface brightness will be ^ 20 mag 
arcsec"^ ((these numbers are nearly independent on whether 
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Figure 13. The final gas surface density profile (tliick solid line) 
of model Lmd2cl2g (see Table 1) is plotted against the initial 
purely exponential profile (thin solid line). The results for the 
isothermal run are shown. Exponential fits to the inner (short- 
dashed line) and outer (long-dashed line) part of the profile are 
shown, with scale lengths of, respectively, 800 pc and 9 kpc (while 
the initial disk scale length was ~ Skpc). The surface density is 
expressed in the units of the simulation 



we consider the case of the models with only a stellar compo- 
nent or we compute surface densities from the models with 
gas assuming that this has been turned into stars). Such 
systems would resemble the red LSBs studied by Bojcrsber- 
gcn, de Blok & Van dcr Hulst (1999). The latter have a low 
surface brightness disk with central red bulges similar in 
structure to those of HSB galaxies but slightly fainter than 
them. Their bulgo-to-disk ratios {B/D) in the I band (which 
should provide a good measure of the actual mass ratios) 
vary between loss than 0.05 to as much as 0.5 and their to- 
tal (disk -I- bulge) B band luminosity is 17 < Mb < 21. The 
final states of our models would have comparable luminosi- 
ties, 18.5 < Mb < 19.5 (for (M/Lb), = 2) and the bulge 
contributes from 15 to 50% of the total mass of the baryons 
(we measure the mass of the bulge as the mass contained 
within the radius at which the slope of the stellar profile 
steepens significantly). A similar evolutionary scenario has 
been recently proposed by Noguchi (2002) to explain giant 
LSB galaxies which bear a resemblance to the red LSBs but 
have a much larger size, bigger B/D ratios and even red- 
der colors (Sprayberry et al. 1995). However, while Noguchi 
starts from a model of a normal HSB galaxy (without gas) 
here we have shown that even galaxies that start with fairly 
low surface density disks can become bar unstable and form 
a bulge provided that their disk is sufflciently massive with 
respect to the halo. The bulges that form in our models 
have both profiles and scale lengths that match those of the 
red, early-type LSB galaxies in Beijersbergen et al. (1999). 
; they are fit by exponential laws and their scale lengths are 
typically ~ 0.1-0.2Rh (Figure 13, 14, and 15). When trans- 
lated into physical sizes, bulge scale lengths are as large as 



0.4-0.9 kpc, i.c significantly larger than the bulges of "nor- 
mal" HSB galaxies of similar luminosity (de Jong & Van 
der Kruit 1994). Being the product of secular disk evolu- 
tion, the bulge scale lengths are correlated with those of 
the surrounding LSB disks, matching another feature of the 
galaxies in Beijersbergen et al. (1999) (for HSB gala:xies a 
similar correlation is indeed found, see de Jong & Van der 
Kruit 1995 and MacArthur, Courteau & Holtzman 2003). 
Also, as already mentioned in section 3.3, disks remain ex- 
ponential but increase their scale length by at least a factor 
of 2, reaching as much as 9 kpc (e.g. Figure 13). Red early- 
type LSBs also have disks with huge scale lengths, bigger 
than their blue counterparts for the same luminosity. Some 
of the red LSBs also exhibit significant spiral structure; this 
appears to be triggered by the bar in our simulations and 
persists even after this has turned into a bulge, especially 
for purely gaseous disks (Figure 7 and 10). 

In our scenario bulges form after a (gaseous) disk has 
already settled into the halo as a result of secular bar evo- 
lution. The bulge would be a younger than the disk in a 
dynamical sense but after a few Gyr it would look redder 
than the latter due to the different mean age of its stars; 
indeed, as we showed above, the different density structure 
of the two components suggests that, while an early burst 
of star formation is plausible for the bulge, the disk would 
undergo a weaker but more prolonged star formation. In 
particular, the spiral arms might excited by the central bar 
would play an important role in the star formation history 
of the disk. Bars and rings are actually present in some of 
the red LSBs, in which case they also have bluer colors in 
their central part; these systems resemble the intermediate 
stages seen in our simulations, while systems with distinct 
bulges arc likely in the late evolutionary stage, and indeed 
they also look redder in the center as if the stellar popula- 
tion there has already undergone significant fading. Future 
detailed observations of color and age gradients throughout 
the disk and the bulge should enable to better judge how re- 
alistic is the scenario proposed here. Simulations including 
star formation and explicit cooling and heating are better 
suited to explore further this issue and will be the subject 
of a forthcoming paper. 

How large a population of LSB galaxies with 
bars/bulges should we expect based on our model? The 
naive expectation is that such galaxies should be rare as 
they originate from systems with disk mass fractions close 
to the upper limit for the baryon fraction (Jaffe et al. 2001). 
On the other end, several authors, by modeling galaxy rota- 
tion curves, find a positive correlation between the spin pa- 
rameter and the disk mass fraction (e.g. Burkert (2000) , Van 
den Bosch et al. (2001), Jimenez, Verde & Oh (2003)), which 
would imply that massive disks are more common among 
LSB galaxies. This correlation is not well understood. Burk- 
ert (2000) has proposed that it might result from a correla- 
tion between disk specific angular momentum (which might 
be different from that of the parent dark matter halo) and 
the disk mass fraction - however at the moment it cannot 
be excluded that the correlation is produced by systematic 
errors in the procedure used to determine the dark halo pa- 
rameters, especially when the rotation curves do not extend 
far enough in radius (Verde et al. 2002). Nonetheless, the 
recent analysys of 400 rotation curves reported in Jimenez 
et al. (2003) indicates that more than 50% of high spin sys- 
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Figure 14. The final gas surface density profile (thick solid line) 
of model Lmdlc4g (see Table 1) is plotted against the intial purely 
exponential profile (thin solid line). The results for the isother- 
mal run are shown. Exponential fits to the inner (short-dashed 
line) and outer (long-dashed line) part of the profile are shown, 
with scale lengths of, respectively, 400 pc and 7 kpc (while the 
initial disk scale length was about 5 kpc). The surfa<;e density is 
expressed in the units of the simulation 



terns (i.e. with A > 0.065, the minimum value considered 
in this paper) have disk mass fractions near 0.1 (L. Verde, 
private communication). These results leave open the in- 
teresting possibility that many LSB galaxies are massive 
enough to become bar unstable and undergo the morpho- 
logical transformation described in this paper. 

Relating the final state of bar unstable but light (fd < 
0.1) gaseous disks to known galaxies turns out to be loss ob- 
vious. The bulge-like component that appears late in these 
systems would have an unusually low optical surface bright- 
ness as it forms out of gas having very low surface densi- 
ties; assuming (M/Lb)^ ~ 4, as typical of old spheroid stel- 
lar populations, these bulges would have a B band surface 
brightness ~ 24 mag arcsec"^. Such a low surface brightness 
is significantly lower than any of those of the objects studied 
in Bcijcrsborgen ot al. (1999). However, although the bulge 
will be fainter than the disk and might be hardly recogniz- 
able in optical bands, it should still stand out when observed 
with longer wavelengths because of its intrinsically higher 
stellar density (see Figure 14). Many blue LSB galaxies or 
LSB dwarfs do not have significant spiral/non-axisymmctric 
structure, appearing rather amorphous or irregular (de Blok, 
Van der Hulst & Bothun 1995); on the basis of our results, 
the simplest interpretation of these systems is that they are 
stable because they have fairly light disks, fd < 0.1 and/or 
halos as highly concentrated as expected in a LCDM model. 
Most of the galaxies whose measured rotation curves are 
nourishing the debate on dark matter cores would fall into 
the latter category; therefore one would be tempted to con- 



clude that bars cannot occur in these systems and hence 
cannot be invoked to explain the present-day structure of 
their inner halo. However, we cannot exclude that the avail- 
able observations lack sufficient resolution and are still hid- 
ing some important clues to the past dynamical histories 
of these objects (see de Blok, Walter & Bell 1999). New 
high-resolution observations of individual blue LSB galaxies 
at long wavelengths will show whether dim, red bulges are 
present in galaxies previously classified as blue, late-type 
LSBs and hence whether bar formation could have taken 
place in the past. Only if such spheroidal components will 
turn out to be extremely rare we will conclude that bar 
formation, and thus bar/halo interactions, can bo neglected 
for these galaxies. Interestingly, the recent observations in J 
and K bands by Galaz et al. (2002) suggest that at least a 
fraction of the systems that appear featureless in the opti- 
cal hide a central red bulge or nucleus when observed in the 
near-IR. Of course, high-resolution images might also reveal 
bar-like distortions. Indeed, some bars are observed in a few 
LSB galaxies included in a recent sample used by Swaters 
et al. (2003) to measure Ha rotation curves; these authors 
claim that the barred systems have slightly shallower halo 
inner slopes compared to the non-barred ones, although they 
admit that non-circular motions due to the bars themselves 
make the determination of the actual rotational velocity very 
hard and questionable in these cases. 

However, even though evidence of present or past bars 
will be found in a number of LSB galaxies, our simulations 
do not provide strong support for a scenario in which bars 
can significantly affect the structure of the dark halo. In fact, 
the bars arising in our simulated galaxies axe rather short, 
from 3 to 8 times smaller than the halo scale radius, (see 
section 2), hence typically shorter than the long massive bars 
assumed in Weinberg & Katz (2002) (but note that the bar 
in model Lmd2cl2b, which is the longest due to the large 
value of the spin i^aramctcr, approaches the size of the bar 
considered in that paper); recent calculations by Sellwood 
(2002) show that, no matter how effective is the transfer 
of angular momentum between the bar and the halo, bars 
up to 4 times smaller than rs do not store enough angu- 
lar momentum to produce the observed large dark matter 
cores. A particularly short bar arises in light disks (model 
Lmdlc4g), its size being > 10 times smaller than r^, hence 
smaller than any of the bars studied by Sellwood (2002). 
Finally, the fact that all the bars in our simulations, and 
especially the gaseous ones, seem to undergo a fast mor- 
phological evolution, also goes in the direction of a modest 
dynamical impact on the dark halo. 



5 SUMMARY 

In this paper we have shown that bar formation in LSB 
galaxies is possible, yet the initial conditions required for 
this to happen apply to only a narrow region of the param- 
eter space made available by currently favored hierarchical 
models of structure formation. In addition, the characteris- 
tics of the bars that eventually form in such systems as well 
as their rapid morphological evolution hardly support the 
idea that bar-halo interactions play a crucial role in shap- 
ing the inner structure of the dark matter halo. Our main 
findings are summarized as follows; 
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• The halo/disk mass ratio within the region where the 
disk lies determines whether the studied LSB models are 
stable or not. This ratio is fixed once both the halo con- 
centration and the disk mass are fixed and sets the degree 
of swing amplification of m = 2 modes. In particular, LSB 
galaxies with disks as massive as 10% of the halo mass can 
become bar unstable for quite a range of concentration val- 
ues. 

• LSB galaxies with (typical) light disks {fd < 0.1) can 
become bar unstable if their halo concentration is as low as 
c = 4 and their disks are essentially gaseous and cold. Such 
low concentrations are rare among LCDM halos. 

• Bars forming in LSB galaxies are usually much shorter 
than the halo scale radius - they are not expected to contain 
enough angular momentum to turn the inner cusp of dark 
matter halos into a core. 

• Both gaseous and stellar bars evolve into central bulge- 
like structures after a few Gyr. The bar and the bulge have 
surface densities up to 2 order of magnitudes higher than 
the surrounding low surface density disks. Radically different 
star formation histories axe thus expected in the central part 
of the galaxy as opposed to the extended disk. 

• The structural properties of the final states (after up to 
10 Gyr of evolution) of bar unstable, massive LSB disks are 
consistent with those of the red, early- type LSB galaxies ob- 
served by Beijersbergen et al. (1999). Our secular evolution 
scenario naturally explains observed correlations like that 
between the disk and bulge scale lengths. 

• Blue LSB galaxies included in samples used to measure 
rotation curves usually appear featureless in the optical; our 
findings imply that that a dim, red bulge-like component 
must be present at their center if they ever formed a bar. 
Future observations of these galaxies in the near infrared 
should reveal whether these hidden bulges exist in most 
LSBs; first attempts in this direction (Galaz et al. 2002) 
suggest that this might be the case. 
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